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We introduce a simple model consisting of two or three coupled one-dimensional channels of par- 
ticles with Yukawa interactions. For the two channel system, when an external drive is applied only 
to the top or primary channel, we find a transition from locked flow where particles in both channels 
move together to decoupled flow where the particles in the secondary or undriven channel move at a 
slower velocity than the particles in the primary or driven channel. Pronounced commensurability 
effects in the decoupling transition occur when the ratio of the number of particles in the top and 
bottom channels is varied, and the coupling of the two channels is enhanced when this ratio is an 
integer or a rational fraction. Near the commensurate fillings, we find additional features in the 
velocity-force curves caused by the slipping of individual vacancies or incommensurations in the 
secondary channels. For three coupled channels, when only the top channel is driven we find a 
remarkably rich variety of distinct dynamic phases, including multiple decoupling and recoupling 
transitions. These transitions produce pronounced signatures in the velocity response of each chan- 
nel. We also find regimes where a negative drag effect can be induced in one of the non-driven 
channels. The particles in this channel move in the opposite direction from the particles in the 
driven channel due to the mixing of the two different periodic frequencies produced by the discrete 
motion of the particles in the two other channels. In the two channel system, we also demonstrate a 
ratchet effect for the particles in the secondary channel when an asymmetric drive is applied to the 
primary channel. This ratchet effect is similar to that observed in superconducting vortex systems 
when there is a coupling between two different species of vortices. 
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I. INTRODUCTION 

There are many systems composed of repulsively in- 
teracting particles with one dimensional (ID) or quasi- 
ID motion, including colloids in narrow channels [l|-[lOl|. 
Wigner crystal states in wires [lll - [l9| and constrictions 
[20] , dusty plasmas in grooves [2l[ , macroscopic charged 
ball bearings in channels [l^l, and vortices in type-II su- 
perconductors confined within narrow strips or channels 
[231427]. In many of these systems, interesting structural 
transitions from ID lines of particles to zig-zag or buck- 
led states can occur d, [13, [l3 ■ There can also be higher 
order transitions from 2 rows to 3 rows of particles or 
transitions to disordered states jS, 23, 24, 27]. Under an 
external drive, this type of system also exhibits a vari- 
ety of dynamical behavior such as ordered or disordered 
motion through constrictions 

0,113,123 or dynamic com- 
mensurability effects [13] . 

Here we propose a simple system consisting of parti- 
cles in two or three coupled ID channels. The particles in 
each channel interact with the other particles in the same 
channel as well as with particles in adjacent channels via 
a Yukawa potential. An external drive applied to only 
one channel produces drag effects on the particles in the 
undriven channels, causing them to move. Our system is 
illustrated in Fig. [T] The two channel system is similar 
to the transformer geometry studied for vortices in two 
superconducting layers where an external drive is applied 
to one (primary) layer and the resp onse of the nondriven 
(secondary) layer is measured [28l - [3C| . If the vortices in 



the two layers are fully coupled, the response of the sec- 
ondary layer is exactly the same as the response of the 
primary layer. If the vortices are only partially coupled, 
the response in the secondary layer is smaller than that 
of the primary layer. The transformer geometry has also 
been studied for vortex systems with multiple layers, such 
as vortices in the strongly layered high-temperature su- 
perconductors [3l| . Drag effects have also been predicted 
for two coupled ID wires containing classical ID Wigner 
crystals when only one of the wires is driven fl^. In 
this case the interaction between particles in neighbor- 
ing wires is repulsive, unlike the attractive interaction 
between vortices in neighboring layers. For the coupled 
ID Wigner crystals there is also a transition from a com- 
pletely locked state, where the response in both wires is 
identical, to a partially locked state, where the response 
of the secondary wire is reduced. Drag effects in coupled 
wire experiments have been interpreted as arising from 
the formation of Wigner crystal states [s^ . 

In this work we consider the effect of changing the ratio 
of the number of particles in each channel on the locking 
or coupling between the channels, with a particular focus 
on ratios that are integers or rational fractions. Com- 
mensuration effects [33] occur when the spacing between 
particles in one channel is a simple rational fraction of 
the spacing between particles in another channel, while 
incommensurations such as vacancies or interstitials ap- 
pear when the two spacings are incommensurable. Com- 
mensuration effects have been studied extensively for sys- 
tems in which a varied number of particles interacts with 
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FIG. 1: Image of the sample geometry. The locations of the 
three channels are indicated by dotted lines. Black dots are 
particles within the channels. The arrow denotes the driving 
force which is applied only to particles within the top channel, 
termed the primary channel p. The bottom undriven channels 
are the secondary channels si and 52- The ratio of the number 
of particles in each channel is -Rsi.p = Nsi/Np and Rs2,p ~ 
Ns2/Np, where Nai and A'^s2 are the number of particles in the 
secondary channels and Np = 16 is the number of particles 
in the primary channel. The spacing between channels d is 
marked in panel (a), (a) Two channels with Rsi,p ~ 1.0 
(b) Two channels with Rsi,p = 0.5. (c) Three channels with 
Rsi,p ~ 1.0 and Ra2,p = 1.0. (d) Three channels with Rai,p = 
1.0 and Rs2,p = 1.5. 



a rigid periodic substrate, such as atoms and molecules 
34| . vo rtices i n superconductors with peri- 



on surfaces 

odic pinning arrays |35l - l37| , and colloids interacting with 
ID [ssl or two-dimensional (2D) optical trap arrays |39|, 
all of which can be viewed as physical realizations of the 
Frenkel-Kontrova model. These studies find that the cou- 
pling to the substrate or the effective pinning of the parti- 
cles by the substrate is strongly enhanced when the ratio 
of the number of particles to the number of substrate 
minima is an integer or a rational fraction, as indicated 
by the appearance of peaks in the critical depinning force 
or enhanced ordering of the particles at the commensu- 
rate fillings. In our system, for the two channel geom- 
etry illustrated in Fig. [T] the particles in the secondary 
channel can be regarded as a distortable or moveable 
periodic pinning substrate for the particles in the pri- 
mary channel, suggesting that enhanced drag or coupling 
could occur when the ratio of the number of particles in 
each channel is an integer or a rational fraction. The de- 
formability of the substrate makes our proposed model 
distinct from Frenkel-Kontrova systems. Additionally, 



driven ID and 2D commensurate-incommensurate sys- 
tems often exhibit numerous dynamic behaviors within 
the incommensurate regimes, such as when localized va- 
cancies or interstitials form soliton-like excitations which 
move more easily than the particles over the substrate 
pol . This suggests that similar phases may be possi- 
ble in the coupled channel drag system we propose here, 
and we show that such phases do appear. We also show 
that when we make the system more complex by adding 
a third channel, a remarkable variety of commensuration 
effects and dynamic regimes occur such as multiple de- 
coupling, recoupling, and slip transitions, all of which 
produce pronounced changes in the velocity response. It 
is even possible to realize negative drag effects where the 
particles in one of the channels move in the direction op- 
posite to that of the applied drive. 

The coupled channels system we propose could be re- 
alized in colloidal systems. The number of colloids in 
the different channels can be controlled readily by opti- 
cal manipulation and the colloids in one channel could 
be driven with an external field, optically, or using mi- 
crofiuidics. Another possible realization of this system 
is in nanowires where ID Wigner crystallization of the 
electrons has occurred; in this case, by altering the elec- 
tron density, the particle lattice spacing in one wire could 
be varied with respect to that in an adjacent wire. Re- 
alizing such a system could have important implications 
for the study of ID Wigner crystals since the appearance 
of commensuration effects would be strong evidence that 
Wigner crystal states are forming. In superconducting 
systems, the density of magnetic vortices is fixed by the 
externally applied magnetic field, so it would be difficult 
to create ID channels that contain different linear den- 
sities of vortices; however, in certain layered systems an 
additional transverse magnetic field can be applied to cre- 
ate a second Josephson vortex lattice which can interact 
with the pancake vortices in the planes [42l - |4^ . It has 
already been shown that using this technique it is pos- 
sible to drive only one of the vortex species and induce 
a drag on the other vortex species [ij, It should 

be possible to study fractional commensurate states in 
such a vortex system by examining how the drag effect 
changes when the ratio of the number of one type of vor- 
tices to the other is varied. A realization of three or more 
channels with varied numbers of particles in each chan- 
nel should again be possible using colloidal systems or 
metallic wires. Further, a superconducting or nanowire 
system could be used in which each layer or channel has 
the same number of particles but differing amounts of 
quenched disorder. We note that there are previous stud- 
ies of colloidal particles in 2D bilayers [45[ where the 
particles in the layers are driven in opposite directions; 
however, these studies focused on an oscillatory order- 
disorder transition, not on the effects of commensuration 
on decoupling or the dynamic phases that we consider 
here for the case of ID coupled channels. 

The paper is organized as follows: In Section II, we 
describe our simulation method and sample geometry. 
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We consider two channels of particles in Section III and 
illustrate a drive-induced decoupling transition for com- 
mensurate channels in Section III A. In Section III B 
we describe the two step decoupling transition that oc- 
curs for incommensurate channels which contain vacancy 
or interstitial sites that can act like a second species of 
particle. The effects of finite temperature and finite size 
appear in Section III C. Section III D shows that the non- 
linear response of the system can be exploited to create 
a ratchet effect, where ac motion in the driven channel 
induces dc transport in the drag channel. In Section IV 
we turn to samples with three channels. We show in 
Section IV A that when the driven channel is commen- 
surate with the neighboring drag channel, four different 
types of coupled and decoupled fiow can occur as the oc- 
cupancy of the second drag channel is varied, including 
regimes of intermittent coupling. In Section IV B, the 
driven channel is incommensurate with the neighboring 
drag channel and we find a complex series of coupling- 
decoupling transitions that produce a significant amount 
of structure in the velocity-force curves. In Section IV 
C, we consider in detail the negative drag that can occur 
at incommensurate fillings when the particles in one of 
the drag channels move in the direction opposite to the 
particles in the driven channel. It is also possible for the 
outer channels to remain coupled while the central chan- 
nel is decoupled, as described in Section IV D. In Section 
IV E we summarize all five of the dynamical phases and 
the negative drag by showing that they can be achieved 
in a single system. The paper concludes in Section V 
with a discussion and summary. 



force is Ff^ = X^j^-i The potential has a 

Yukawa or screened Coulomb form of 

ViR.,) = |^Le-««- (2) 

with Eq ~ Z*^/4:TTe£oao, where e is the solvent dielec- 
tric constant, Z* is the effective charge, and 1/k is the 
screening length. For colloidal systems, the length scale 
ao is on the order of a micron. We measure forces in 
units of Fq = Eo/aQ and time in units of r = ti/Eq. In 
a typical case, the distance between particles in channel 
p is a = 2.25 while the distance between adjacent chan- 
nels is d = 1.125a and the screening length is 1/k — 2d, 
which is long enough to ensure strong coupling between 
the particles in all three channels. The driving force 
= Fdx is applied only to all the particles in the 
primary channel. We increase Fu from zero in small in- 
crements of SFd, holding the drive at constant values 
for a fixed time interval during which we measure the 
velocity of the particles in each channel. We have care- 
fully checked that our waiting times are long enough to 
eliminate transient effects. We use SFd = 0.001 and a 
wait time of 10^ simulation time steps. We impose pe- 
riodic boundary conditions in the x-direction along the 
length of the channels. The velocity of the particles in 
the primary channel is given by Vp and the particle ve- 
locity in the secondary channels is given by Vsi and Vs2- 
We normalize all the velocities by the number of particles 
in each channel, Vp = N-^ ^f" v„ Vsi = Ef°' v„ 
and = iV72i Ef"v,. 



II. SIMULATION 

We model interacting Yukawa particles confined to 
move along ID channels as illustrated in Fig. [TJ Each 
particle interacts with other particles in the same channel 
and with particles in adjacent channels. The separation 
between channels is d = 2 and, unless otherwise noted, 
there are Np = 16 particles in the driven or primary 
channel with a lattice spacing a, where L is the length 
of the channel. The particles in the primary channel p 
are coupled to an applied external driving force Fo- For 
a two channel system, the drag or secondary channel si 
contains N^i particles, and the commensurability ratio is 
Rsi,p = Nsi/Np. In a three channel system such as that 
shown in Fig.[IJc,d), the additional secondary channel S2 
is adjacent to si but not to the primary channel p, and it 
contains Ns2 particles, giving a commensurability ratio 

of Rs2,p = Ns2/Np. 

The particle motions evolve under overdamped dynam- 
ics where the colloids obey the following equation of mo- 
tion: 



Here 77 is the damping constant, is the location of 
particle i, and the repulsive particle-particle interaction 



III. TWO CHANNEL SYSTEMS 

A. Coupling-Decoupling Transitions for 
Commensurate Channels 

We first focus on the two channel system at commensu- 
rate filling with a particle ratio of Rsi,p = 1.0 and with 
d/a = 0.67. In Fig. [Ha) we plot Vp and Vsi together 
versus Fo. Both Vp and Vsi increase linearly with Fd 
for low Fd and have identical values, indicating that the 
motion in the two channels is locked. At Fd = 2.125, 
we find a transition to a partially decoupled state where 
Vsi monotonically decreases with increasing Fd while Vp 
continues to increase with Fd at a rate faster than the 
linear increase that occurred below the transition. The 
particles in si are not completely decoupled from the 
primary channel since they still exhibit a nonzero veloc- 
ity; since the particles in si do not experience a driving 
force, they can move only due to interactions with the 
particles in p. Just above the decoupling transition, Vp 
increases with a square root form. The general shape 
of the velocity force curves in Fig. [UJa) is the same as 
that of the current- voltage curves obtained for supercon- 
ducting transformer geometries [H-fsH , where the vortex 
velocities are proportional to the voltage and the applied 
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FIG. 2: (a) The average velocity in the primary channel 
Vp and the secondary channel Vsi vs Fd for a two channel 
system with Rsi,p ~ 1.0 and d/a — 0.67. The channels are 
locked for low Fd in the regime where Vp and Vai increase 
linearly. At Fd ~ 2.125, there is a transition to a partially 
decoupled state where the particles in the secondary channel 
begin to slip, producing a decreasing Vsi- (b) The decoupling 
force Fc force vs d/a for the system in (a) with Rai^p — 1, 
fixed d = 2.0, and a altered by changing the particle den- 
sity. The decoupling transition drops to lower values of Fd 
as the particle density increases, (c) Vp (upper curve) and 
Vgi = Np^ (lower curve) vs Fd for the system in (a) 

but with Rai,p = 0.92. Unlike the commensurate case in (a), 
the initial unlocking phase above Fd = 1.04 is associated with 
the slipping of vacancies in si, while at Fd = 2.25 there is a 
second unlocking transition above which all the particles in 
si slip with respect to the particles in p. (d) Vp (upper curve) 
and Vsi (lower curve) vs Fd for Rsi,p ~ 1.08, where the slip- 
ping at low drives is due to the presence of incommensurate 
particles. 



current is proportional to the external force on the vor- 
tices in the primary channel. The current-voltage curves 
in the superconducting transformer system indicate that 
there is a drive-induced decoupling transition of the vor- 
tices in adjacent layers. In the vortex system, the vortex- 
vortex interaction between layers is attractive, so it is 
more intuitive why a finite vortex mobility persists in 
the secondary channel above decoupling. The results in 
Fig. [2ja) indicate that even when the particle-particle 
interactions are purely repulsive, a finite velocity in the 
secondary channel can be maintained at drives above the 
decoupling transition. 

We measure the decoupfing or unlocking force as a 
function of d/a in a two channel system with Rsi.p = 1.0. 
The result is plotted in Fig. ^h), where we fix d = 2.0 
and vary a by changing Np and A^si. Here Fc decreases 
with increasing particle density. A similar effect occurs in 
layered vortex systems, where for higher fields or higher 
vortex densities the coupling between the layers is grad- 
ually reduced fsij . For the Yukawa system this effect can 



be attributed to the reduced size of the periodic potential 
that the particles in si experience from the particles in 
p. Once the primary channel is in motion, the particles 
in si shift to positions that are slightly behind the driven 
particles. As the drive increases, the size of this shift in- 
creases until the particles in p slip more than 0.5a ahead 
of the particles in si, producing the partial decoupling. 
When the particle density increases, the amount of shift 
required to pass the position 0.5a decreases since a de- 
creases with increasing particle density in both channels. 



B. Dynamics and Commensurability 

In Fig. ^c) we plot Vp and V^i for a two channel sys- 
tem with i?si,p = 0.92, where the number of particles 
in si is smaller than the number of particles in p. Here 
V^i — Sf^^^ velocity of the particles in si 

normalized by Np, the number of particles in p. For low 
Fd the channels are locked and all of the particles in the 
system move at the same velocity. The slope of V^^ versus 
Fd is slightly smaller than the slope of Vp versus Fd in 
Fig.H^c) due to the fact that Nsi < Np. At Fd = 1.04 we 
observe a transition to a partially coupled state; however, 
this transition occurs at a drive well below the decoupling 
transition Fc — 2.125 shown for the commensurate sys- 
tem in Fig. [2Ia). Additionally, just above Fd = 1.04, 
Fig- IHc) indicates that the velocity-force curve does not 
have the characteristic square root shape found close to 
Fc in Fig. [^a). For Fd > 1.04, V^i continues to increase 
with increasing Fd but with a smaller slope than in the 
locked regime. A second decoupling transition appears at 
Fd = 2.25. For Fd > 2.25, V^^ decreases with increas- 
ing Fd and there is also a corresponding increase in the 
slope of Vp. The second decouphng transition occurs at 
a drive close to the value Fc = 2.125 where decoupling of 
the commensurate system occurs, as shown in Fig.^^a). 
This indicates that the second decoupling transition for 
the incommensurate system is the same as the sole de- 
coupling transition found in the commensurate system, 
where all the particles in si begin to slip with respect to 
the particles in p. The two step decoupling transition for 
the incommensurate system appears due to the presence 
of vacancies in si. At commensuration, all of the par- 
ticles in si are located within potential minima created 
by the spacing of the particles in p. Below commensura- 
tion, a fraction of the sites in this periodic potential are 
empty, producing effective vacancies in si. In the locked 
phase at Fd < 1.04, all the particles in si move at the 
same velocity as the particles in p. At the first decoupling 
transition, the vacancies in si begin to slip with respect 
to the particles in p. This can be viewed as a dcpinning 
transition. Every time a vacancy slips, only one of the 
particles in si slips with respect to p while the remaining 
particles in si stay locked with p. As a result, most of 
the particles in si continue to increase in velocity with 
increasing Fd. For drives above the second decoupling 
transition, all of the particles in si slip with respect to p 
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and the slipping is no longer dominated by the motion of 
vacancies. 

In Fig. [2jd) we plot Vp and V^^ versus Fd for the 
same system in Fig. [2][c) but with Rsi,p = 1.08, where 
Nsi > Np so that a few incommensurate particles ap- 
pear in si. The overall shape of the velocity- force curve 
in this case is very similar to that for Rsi.p = 0.92 shown 
in Fig. Hfc), with a first decoupling occurring at a lower 
drive of Fd = 0.6 than that for i?si,p = 0.92, and a sec- 
ond decoupling transition occurring close to F]j — 2.0. 
Here the incommensurations in si form doubly occupied 
sites in the periodic potential created by the particles in 
p. At Rsi^p = 0.92 when there are vacancies in si, slip- 
ping of a particle adjacent to a vacancy occurs because 
the particle is able to move closer to the barrier sepa- 
rating two minima in the periodic potential. This is be- 
cause the force the particle experiences on one side from a 
neighboring particle in si is not compensated due to the 
missing particle at the vacancy site. As a result, there 
is an extra force of the order of FPP{a) on the slipping 
particle, where a is the lattice constant of the particles 
in p. For the doubly occupied sites at Rsi^p — 1.08, a 
similar situation occurs; however, the slipping particle in 
si is located at a doubly occupied site and feels an un- 
compensated force from the other particle located within 
the same site. The extra force in this case is Fpp{a'), 
where a' < a in order for the site to be doubly occu- 
pied. Thus, Fpp{a') < Fpp(a), so the initial decoupling 
transition occurs at a lower value of Fu for samples with 
Rsi,p > 1 that have incommensurations than for samples 
with Rsi^p < 1 that contain vacancies. 

The appearance of multiple decoupling transitions just 
below and above commensuration and only one transition 
at commensuration is similar to the single and multiple 
depinning transitions observed in vortex systems :4Q] and 
colloidal systems j46i] with periodic potentials at and near 
commensuration. In these 2D systems, at the match- 
ing filling of 1 particle per substrate minimum there is a 
single transition from a pinned state to a flowing state, 
while at fillings slightly away from commensuration, well- 
defined vacancies or interstitial particles appear which 
are highly mobile and depin at a lower external drive 
than the commensurate particles. The 2D systems are 
generally more complicated and allow for more than two 
depinning transitions near but not at commensuration 
[40]; however, as in our ID case, it is the presence of two 
types of particles, the commensurate particles and the 
interstitial or vacancy sites, that produce the multiple 
depinning transitions. There are some important differ- 
ences between our two channel system and the 2D vortex 
and colloidal systems. Our system contains no fixed pe- 
riodic substrate so there is no pinned phase; however, 
there is a moving fully coupled state which is analogous 
to the pinned state. The regime in which the vacancies 
or incommensurations slip is then analogous to the de- 
pinning transitions of interstitials or vacancies, and the 
high driving phases at which all the particles are slipping 
corresponds to the completely depinned regime in the 2D 
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FIG. 3: (a) The velocity n of a single particle in si vs time 
for the commensurate system in Fig. [2ja) at Rsi,p = 1.0 and 
Fn — 2.25. (b) The Fourier transform S{f) of the signal 
in panel (a) shows a well-defined characteristic frequency, (c) 
v{t) for a system with Rsi,p ~ 1.08 at Fd = 0.65. (d) The cor- 
responding S{f) shows that there are two frequencies present, 
(e) v{t) for a system with Rsi,p = 0.92 at Fd = 1.1. (f) The 
corresponding S{f) shows the presence of two frequencies. 



systems. 

In order to show more clearly that the particles in si 
are experiencing a periodic potential produced by the 
particles in p and that there are two effective types of 
particles in si away from commensuration, in Fig. [3l^a) 
we plot the time trace of the velocity v{t) of a single 
particle in si for the system in Fig. HJa) at Rsi^p = 1.0 
in the locked phase at Fu = 2.25. The value of v{t) is 
nearly constant except during the periodic slip events, 
during which v drops briefly below zero indicating that 
the particle temporarily moves backwards. In Fig. 
we show the Fourier transform S{f) of v{t) highlighting 
the presence of a single characteristic frequency deter- 
mined by the slipping events. In the locked phase, there 
is no high-frequency oscillation of the velocity of any of 
the particles. In Fig. ^c) we plot v{t) for the system 
with R.i^p = 1.08 from Fig.llfd) at Fd = 0.65 where the 
channels are not completely locked but where V^i is still 
increasing with increasing Fu, and in Fig. [3l^d) we show 
S{f) for the same data. There are now two frequencies 
present. The lower frequency is produced by the same 
slipping events that occurred for the commensurate sys- 
tem in Fig. [3l^a,b), while the higher frequency originates 
from the motion of the incommensuration through si. 
In Fig. [31[e,f) we plot v{t) and S{f) for a sample with 
Rsi,p = 0.92 at Fd — 1-1. Again, there are two char- 
acteristic frequencies. The lower frequency is associated 
with the same slipping events shown previously, while the 
higher frequency is produced by the motion of a vacancy 
through si, rather than by the motion of an incommen- 



FIG. 4: Vp (upper curve) and Vsi (lower curve) vs Fd for 
the two channel system from Fig. [2ja). (a) At Rsi,p ~ 0.5, 
there is a single transition from the completely coupled state 
to the partially decoupled state, (b) At Rai,p = 0.58 there are 
two transitions. The first decoupling transition occurs when 
the incommensurations begin to slip while the other particles 
in si remain locked with p. (c) At Rai^p = 2.0 there is a 
single transition out of the locked phase, (d) At Rai.p = 1.92 
(upper Vsi curve) and Rbi,p = 2.08 (lower Vsi curve) there 
is no locked phase, but there is still a local maximum in Vsi 
which is higher for Rai^p = 1.92. 



suration. 

The two stage decoupling process is most pronounced 
for fillings close to Rsi,p — 1.0, but the same effects 
appear near certain fractional ratios. For example, in 
Fig. Ufa) we plot Vp and Vsi versus Fd for a sample with 
Rsi.p = 0.5 which has a single sharp decoupling transi- 
tion at Fd — 1.9. Just above this filling at Rsi,p = 0.58, 
shown in Fig. Ul^b), there is an initial decoupling transi- 
tion of the incommensurations near Fd = 1.5 into a state 
where Vsi still increases with increasing Fd but with a 
greatly reduced slope. A second decoupling transition 
appears at Fd = 1.8, and above this drive Vsi decreases 
with increasing Fd ■ In this case, incommensurations ap- 
pear with respect to the particle configuration that oc- 
curs at Rsi^p = 0.5. There is a single frequency associ- 
ated with the motion of the particles in Si for i?si_p = 0.5, 
while for the incommensurate case of Rsi,p = 0.58, two 
frequencies are present. This trend persists for higher 
values of Rsi.p as shown in Fig.Hl^c) at Rsi^p = 2.0. Here 
there is a single sharp decoupling transition, while just 
below and just above this filling at = 1.92 and 

Rsi,p = 2.08, Fig. \M.<i) shows that the locking phase is 
absent but that a strong local maximum in V^i appears 
at Fd = 0.19 for Rsi,p = 1.92 and at Fd = 0.125 for 
■Rsi,p = 2.08. 

By performing a series of simulations for varied Rsi,p, 
we determine the location Fc of the transition from com- 
plete locking to a decoupled state and map out where the 
commensuration effects occur. In Fig. [SJa), we plot the 




FIG. 5: The force Fc at the transition from the locked to 
unlocked phase vs Rsi,p in two channel samples with different 
values of a. (a) At d/a — 0.44, commensuration peaks occur 
at Rsi,p = 1.0, 2.0, and 3.0. Fractional peaks and anomalies 
appear at Rsi,p = 0.5, 1.5, and 2.5. (b) At d/a — 0.67, there 
are commensuration peaks at Rai^p = 0.5, 1.0, 1.5, and 2.0. 
(c) At d/a = 1.0, the strongest commensuration peaks appear 
at Rsi.p = 1 and 0.5. 



decoupling force Fc versus Rsi.p for a two channel sample 
with d/a = 0.44, which falls in the strong coupling regime 
in Fig. [2jb). There are peaks in Fc at Rsi,p = 1.0, 2.0, 
and 3.0, along with submatching peaks at Rsi.p — 1/3, 
1/2, and 2/3. Additionally, weaker anomalies appear at 
Rsi.p = 1-5 and 2.5. In Fig. [S]Jb) we show Fc versus 
Rsi.p for a sample with d/a = 0.67. The value of Fc at 
the commensurate filling of Rsi.p = 1.0 is lower for the 
d/a = 0.67 sample than for the d/a = 0.44 sample. Fig- 
ure [5Ub) also has clear peaks in Fc at Rsi,p — 2.0, 1.5, 
and 0.5, while above Rsi.p = 2.0 within our resolution 
there are no peaks or regions where the system is locked. 
In the regions with Fc = where the locked phase is 
absent, the second decoupling transition still appears at 
higher drives and can be detected as the point at which 
Vsi changes from increasing to decreasing with increasing 
Fd- For higher particle densities and fixed d, the com- 
mensurability effects still persist as shown in Fig.[5jc) for 
d/a = 1.0. Here, peaks in Fc occur at Rsi.p = 0.5, 1.0, 
and 2.0. 

The appearance of the commensuration effects at in- 
teger and fractional fillings suggests that this system ex- 
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FIG. 6: Vsi vs Fd for the system in Fig. [5jb) with d/a — 
0.67. (a) Rsi,p = 0.562, 0.625, 0.6875, and 0.75, from top to 
bottom. Inset: Detail of the -R3i,p = 0.6875 curve from the 
main panel, (b) Rsi.p = 0.8125, 0.875, 0.9375, and 1.0, from 
bottom to top. (c) R,i,p = 1.0625, 1.125, 1.1875, and 1.25, 
from top to bottom. 



hibits the same behavior found for the depinning of re- 
pulsively interacting particles on a ID fixed periodic po- 
tential; however, there are several differences between the 
two systems. For particles on a fixed periodic potential, 
the depinning force at fields where the particle-particle 
interactions cancel due to symmetry equals the maximum 
value of the pinning force Fp so that Fc — Fp at fillings 
1/12, 1/8, 1/6, 1/4, 1/2, and 1.0. For the drag system 
shown in Fig. [SJ this does not occur and there is even a 
trend for Fc to increase at the lowest fillings. This is be- 
cause the substrate potential created by the particles in 
p is not fixed but can distort since the particles in either 
channel can shift. At Rsi,p = 1.0, the periodic potential 
is fairly rigid due to the matching of the particle posi- 
tions in p and si, and any distortion of the particles in 
p is energetically unfavorable. In contrast, at very low 
fillings such as Rsi,p = 0.125, the particles in p distort 
near the locations of the particles in si in order to cre- 
ate a localized lowering of the density in p above each 
particle in si. As a result, the particles in si no longer 
experience the same periodic potential from p that was 
present for the commensurate case of Rsi,p — 1.0. Even 
at Rsi^p = 0.5, the particles in p can distort, reducing 
the strength of the coupling to the particles in si. 

In order to better understand the changes in dynamics 
at the different fillings, in Fig. [5] we plot T4i as a function 
of Fd for varied Rsi,p in a system with d/a = 0.67. At 
Rsi.p = 0.5, a single decoupling transition occurs and Vsi 



is a monotonically decreasing function. For = 0.562 
and 0.625, shown in Fig.[S]Ja), there is a clear double peak 
structure in Vsi with one peak falling at the depinning of 
the incommensurations and the second peak appearing at 
the unlocking transition. At i?si,p = 0.6875 in Fig. [6ja), 
there is now a three peak structure in Vsi ■ The first peak, 
shown in the inset of Fig. [ni[a), falls at the transition out 
of the completely locked phase at Fd = 0.11. The second 
and largest peak is at F^ = 0.3, while a third broad peak 
also appears that is centered at Fjj — 1.45. The broad 
peak is the remnant of the second peak in V^i found 
for Rsi^p = 0.562 and 0.625; with increasing i?si,p, this 
peak broadens and the center shifts to higher values of 
Fd- For 0.11 < Fd < 0.3, the particles in si are almost 
completely locked but there is a single incommensuration 
which has begun to slip. For = 0.75, the initial peak 
is lost and the decoupling transition peak now falls at 
Fd — 0.11. There is also a very broad maximum centered 
at Fd = 4.0. Another interesting feature is that at higher 
Fd such as at Fd — 6.0, Vsi for Rsi.p — 0.75 is higher 
than Vsi at the lower values of Rsi,p, even though at low 
Fd Rsi.p showed the lowest value of Vsi. This suggests 
that at high values of Fd, additional drag is produced 
by the interaction between the incommensurations in si 
and the particles in p. 

In Fig.[S]Jb) we plot Vsi versus Fd for Rsi.p = 0.8125, 
0.875, 0.9375, and 1.0. The maximum value of Vsi 
increases as Rsi.p increases toward Rsi.p — 1.0 and 
the broad maximum in Vsi sharpens and shifts toward 
lower Fd- Here, Rsi.p — 1.0 has the lowest value and 
Rsi.p — 0.8125 has the highest value of Vsi at Fd = 6.0. 
In Fig. ^c) we show Vsi versus Fd for Rsi.p = 1.0625, 
1.125, 1.1875, and 1.25. For these values of Rsi.p there 
is no completely locked phase; however, there is still a 
peak feature in Vsi for Rsi.p = 1.0625 and 1.125 which 
broadens and shifts to higher Fd for increasing Rsi.p- 
At Rsi.p — 2.0 a locked phase reappears and the shape 
of Vsi versus Fd is very similar to the curve shown for 
Rsi.p — 1-0. 

In Fig. [Tl^a) we plot Vsi versus Rsi.p at a fixed drive 
of Fd — 6.0 for the system in Fig. [S] with d/a — 0.67. 
Larger values of Vsi indicate that the particles in si are 
exerting a larger drag on the particles in p. Here, peaks 
fall at Rsi.p = 0.45 and Rsi.p = 0.8, rather than at the 
values Rsi.p = 0.5 and Rsi.p = 1.0 where peaks appeared 
in Fc in Fig. ^h)- This shows that at high drives, the 
drag by the si particles is the most effective away from 
the commensurate fillings. The curves shown in Fig. [S] 
indicate that if Vsi were measured at a lower value of Fd, 
the peaks in Fig. [TJa) would shift closer to Rsi,p = 0.5 
and 1.0. 

In vortex systems with 2D periodic pinning arrays, ex- 
periments have shown that the pinning is enhanced at 
the matching fields as indicated by dips in the resistiv- 
ity for low applied drives. For vortices that are strongly 
driven, however, the resistivity dips were found to shift 
away from the integer matching fields [43|. The inter- 
pretation was that at the matching fields in the highly 
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FIG. 7: Vai vs Rsi,p at Fn ~ 6.0 for the two channel system 
with d/a — 0.67. Here the peaks appear not at the commen- 
surate fields of Rsi,p = 1.0 and Rsi.p — 0.5 but at Rsi,p = 0.8 



and Rsi,p = 0.45. (b) K'J 



lsi,p obtained for each filling 



at the Fd where Vai reaches its maximum value. The peak 
centered at -Rsi.p = 1.0 is much broader than the peak in Fc 
at _Rsi,p ~ 1.0 shown in Fig. [IJb). 
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FIG. 8: (a) Fc vs Rsi,p for the system with d/a — 0.67 
from Fig. [5jb) of length L (diamonds), 2L (squares, curve 
shifted up by 0.4), and 4L (circles, curve shifted up by 0.8). 
For the larger systems, higher order fractional peaks appear 
at Rsi,p = 1/4, 1/3, 1/2, 2/3, 3/4, and 3/2. (b) The same 
data plotted without vertical shifts, for system sizes of L (di- 
amonds), 2L (squares), and 4L (plus signs). Connecting lines 
are drawn only for the 4L system. The curves overlap exactly 
and the values of Fc are unaffected by system size. (c,d) Vsi 
vs Fd for the 41/ system from (a) at (c) Rai^p — 1.0 and (d) 
J, = 0.896 for temperatures of T = 0.0, 0.22, 0.88, 2.0, and 
4.5, from top to bottom. The decoupling transition drops to 
lower Fd with increasing T, while the drag effects persist up 
to high temperatures. 



driven system, the vortices form a very ordered moving 
commensurate state, while at the incommensurate fields 
the moving state is not as well ordered and thus the ef- 
fectiveness of the pinning increases away from commen- 
suration at high drives. Although the disorder in the in- 
commensurate state causes the system to begin slipping 
at a lower drive for incommensurate fields, at high drives 
the disordered state experiences more fluctuations than 
the ordered state which induce some additional drag. 

In Fig.[7l^b) we plot V™""^ versus Rsi,p- Here the mea- 
surement of VJ^'^^ is performed not at a fixed Fd but at 
the Fjj where Vsi reaches its maximum for each value of 
Rsi,p- In this case, a strong peak in V^""^ appears at 
Rsi,p = 1-0. This peak is wider than the peak in Fc at 
Rsi,p = 1-0 in Fig. mb) due to the fact that the maxi- 
mum value of Vsi increases as i?si,p = 1-0 is approached, 
as shown inFig.|ni;b,c). 



C. Finite Size and Temperature Effects 

To determine whether further higher order submatch- 
ing effects in Rsi,p can be resolved for larger systems 
and whether the values of Fd at which the unlocking 
transitions occur change with system size, we consider 
the system at d/a = 0.67 from Fig. [SJb) and analyze Fc 
for samples of size 2L and AL. Here we hold a fixed by 



increasing A'^p to 32 and 64, respectively. In Fig. [51Ja) 
we plot Fc versus Rsi,p for samples of size L, 2L, and 
4L, with the curves shifted vertically for clarity. In the 
larger samples, there are clearly fractional peaks falling 
at Rsi^p = 1/4, 1/3, 1/2, 2/3, 3/4, and 3/2. The 4L sam- 
ple even shows some evidence of a peak at Rsi.p = 1/8. 
We expect that for even larger systems, even more frac- 
tional peaks will appear but that the higher order peaks 
will be increasingly weak in size, similar to the behavior 
of fractional peaks observed in other systems such as vor- 
tices on periodic substrates [3^[36l]. In Fig. |Sl^b) we plot 
the same data without vertical shifts to show that the de- 
pinning thresholds for the three systems overlap exactly; 
only the resolution is changed by the system size. We 
find no changes in the velocity-force curves as the size 
of the sample is increased, indicating that the system 
sizes we are studying capture the essential behavior. We 
also find a similar lack of dependence on sample size for 
the three layer systems that are described in Section IV. 
We note that for commensurate- incommensurate systems 
such as the Frenkel-Kontrova model [33], submatching 
effects theoretically occur for all rational values oi m/n, 
where m is the number of particles and n is the number 
of substrate minima. In the Frenkel-Kontrova model, 
true incommensurate behavior occurs only for systems 
of infinite size at irrational filling ratios. In our system 
the higher order submatching effects are destroyed due 
to the fact that we do not have a fixed substrate; in- 
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stead, the effective substrate experienced by the parti- 
cles in one channel due to the presence of particles in 
the neighboring channel is able to distort. As a result, 
our system does not map directly onto commensurate- 
incommensurate systems such as the Frenkel-Kontrova 
model, although it displays several similarities with such 
models as we have shown. 

In an experimental realization of the system we pro- 
pose, such as with colloids confined to channels, ther- 
mal effects will be present. To test the stability of the 
different regimes we observe against thermal perturba- 
tions, we have performed simulations with the AL sys- 
tem with d/a = 0.67 at i?si,p = 1-0 and Rsi,p = 0.896 
at finite temperature. We use the same procedure em- 
ployed in previous works to model the thermal fluctu- 
ations [il]. We add a Langevin noise term to the 
equation of motion with the properties {Ff{t)) = and 
{Ff{t)Fj'{t')) = 2rikBT5,j5{t-t'). Fig. IH^c) illustrates 
Vsi versus Fd at Rsi.p — 1.0 and Fig. [8jd) shows V^^i 
versus Fd at Rsi^p = 0.896 for T = 0.0, 0.22, 0.88, 2.0, 
and 4.5. As the temperature increases, the value of Fc 
decreases until for T > 2.0 the locked phase has almost 
completely vanished; however, drag effects on the sec- 
ondary channel continue to persist up to much higher 
temperatures. For Rsi,p = 0.896, the locking phase is lost 
at lower T than for Rsi^p — 1.0 due to the fact that the 
effective incommensurations are more mobile and hence 
require a smaller level of thermal fiuctuations to escape 
from the potential minima. These results show that the 
drag and locking features described for the zero tempera- 
ture system should persist under finite temperature pro- 
vided that the thermal fiuctuations are not excessively 
strong. 

Possible experimental realizations of the two channel 
system include modified versions of the colloidal exper- 
iments which have already been performed on coupled 
one-dimensional channels Q • Colloidal systems are sub- 
ject to thermal fiuctuations and hydrodynamic interac- 
tions which can arise in the surrounding fiuid. We showed 
above that the dynamic phases are robust against moder- 
ate thermal fluctuations. There is ongoing discussion re- 
garding how the inclusion of hydrodynamic effects would 
impact the dynamics of driven colloidal systems. Re- 
cent two-dimensional simulations of an electrophoreti- 
cally driven charged colloidal system similar to the one 
we consider showed that when the charge on the colloids 
are sufficiently strong, the dynamical behavior of the 
system is not altered by the addition of hydrodynamic 
interactions [1^. Due to the good agreement that has 
been found between numerous simulations of driven col- 
loid systems in which hydrodynamic effects are neglected 
and the actual behavior of driven colloids in experiment, 
we expect that at least some of the features that we de- 
scribe should be observable in a colloidal realization of 
this system. Another possible realization would be to 
generalize the recent experiments performed with dusty 
plasmas interacting with a one-dimensional groove to cre- 
ate what are termed Yukawa chains '211. If more than 
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FIG. 9: Inset: Schematic of the ac force applied to p. Each 
period r is divided into two parts. A force Fa is applied in the 
positive direction for a duration ta, then a force Fb is applied 
in the negative direction for a duration tb, with the condition 
that Fb/Fa = ta/tb = 4.0. Main panel: Induced dc velocity 
Vai averaged over multiple ac drive periods vs Fa for a two 
channel system with d/a = 0.67 under the applied ac drive 
shown schematically in the inset. Upper curve: Rai,p = 1.0; 
middle curve: Rsi,p = 0.75; lower curve: Rsi,p = 1.25. Here 
the ratchet effect reaches its max;imum value for Rs\,p ~ 1.0. 
Within this range of Fa, Vp — when averaged over an ac 
drive period. 



one groove were created in the substrate, it should be 
possible to couple two or more of the Yukawa chains, to 
use a laser to drive one of the chains, and then to analyze 
the response of the secondary chain. In this case inertial 
effects could modify the behavior since dusty plasma sys- 
tems are generally not in the overdamped limit. 



D. Ratchet Effect With ac Drives 

We next show that when the particles in p are driven 
with an ac drive, it is possible to generate a net dc motion 
of the particles in si or a ratchet effect. Ratchet effects 
produced by applied ac drives have been studied exten- 
sively in systems of particles interacting with asymmet- 
ric substrates [13]; however, it is also possible to create a 
ratchet effect in the absence of an asymmetric substrate 
when the ac drive has certain asymmetries and when the 
response of the system is nonlinear [i^ - H -lslj. This type 
of ratchet has been realized in systems with two inter- 
acting species of superconducting vortices such as when 
Josephson vortices couple to pancake vortices j43 - |4^ , as 
well as in interacting binary colloidal systems where only 
one colloid species couples to an external driving field 
and produces a rectification of the other colloid species 

Here we consider an ac square drive applied only to 
p. The period r of the square drive is divided unevenly 
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into two parts as illustrated in the inset of Fig. |9l In 
part A, we apply a force = i^^x in the positive di- 
rection for a duration ta, while in part B we apply a 
force Fd = —Fb^ in the negative direction for a dura- 
tion tb — T — TA- In selecting Fa and Fb, we impose 
the condition Fata — Fbtb = so that there is no net 
dc drive. If the response of the system is perfectly lin- 
ear, this drive will not generate a net dc motion of the 
particles in either channel. On the other hand, if the 
coupling between si and p is nonlinear, it is possible to 
induce a dc motion of the particles in si by applying this 
ac drive to the particles in p. If both Fa and Fb are 
below the first decoupling transition Fc, the motion of 
the particles in both channels is completely locked, the 
response is perfectly linear, and there is no ratchet effect. 
If Fa < Fc and Fb > Fc, a net dc drift of the particles in 
si will occur since the particles in si remain completely 
locked with the particles in p during part A of the drive 
cycle, but during part B of the cycle the particles in si 
are partially decoupled and do not move all the way back 
to their starting position by the end of the cycle. 

In Fig. [5] we illustrate the ratchet effect which produces 
a finite positive value of Vsi under the ac drive described 
above. We fix Fb/Fa = 4.0 and d/a = 0.67, and plot the 
time-averaged l^i versus Fa for Rsi,p = 0.75, 1.0, and 
1.25. For low Fa, Vsi starts small but rapidly grows with 
increasing Fa, reaching a sharp peak for Rsi^p = 0.75 and 
Rsi,p — 1-0. As Fa increases above this peak, Vsi grad- 
ually decreases with increasing Fa since the drag effect 
becomes smaller for higher drives as shown in Fig. [D^a) . 
For Rsi,p = 1.25, the ratchet effect is strongly reduced 
but still persists, indicating that the ratchet effect should 
be a robust feature for all fillings. In all cases there is no 
induced dc flow of the particles in p. Our system can be 
regarded as containing two species of particles: the di- 
rectly driven particles in p, and the undriven particles in 
si that experience a drag from the particles in p. It would 
be very interesting to look for a similar ratchet effect in 
coupled quantum wires in the regime where Wigner crys- 
tallization may be occurring. This could be achieved by 
applying an ac drive of the type illustrated in the inset of 
Fig. inito one wire and determining whether a dc response 
is induced in the second wire. 
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FIG. 10: The velocities Vp, Vsi, and Vs2 vs Fd for a three 
channel system with d/a = 0.67, Rsi,p ~ 1.0, and varied 
Rs2,si- (a) At Rs2,si = 1.0 there is a single transition from 
the locked region I to region II where the particles in si and 
S2 remain locked with each other but partially decouple from 
the particles in p. (b) At i?s2,3i = 1.16 there is a transition 
from the locked region I to region III where the particles in si 
and p lock together but the particles in S2 partially decouple. 
In region IV all three channels are unlocked, while at high 
Fd the system enters region II when the particles in si lock 
with the particles in S2 but are partially decoupled from the 
particles in p. (c) At Rs2,3i = 1-5 the transition between 
regions IV and II occurs at a much higher value of Fd ■ 



IV. THREE CHANNEL SYSTEMS 

A. Coupling-Decoupling Transitions for Partial 
Commensuration 

We next consider a system with three channels of par- 
ticles where only the top channel is subjected to a driving 
force. We measure the velocities in each channel, denoted 
by Vp, Vsi, and Vs2, for particle ratios of i?si,p — Nsi/Np, 
Rs2,p = Ns2/Np, and i?s2,si = Ns2/Nsi. For the com- 
mensurate case when all channels contain the same num- 
ber of particles, Rsi,p = Rs2,p = Rs2,si = 1-0, the behav- 
ior is the same as in the two channel case at commensura- 



tion. There is a single decoupling transition from region 
I, the completely locked phase, to region II, where the 
particles in si and S2 remain locked with each other but 
are partially decoupled from the particles in p. This is 
illustrated in the plot of Vp, Vsi, and Vs2 versus Fo in 
Fig.fTUTa') for a sample with d/a = 0.67, Rsi,p = 1.0, and 
Rs2.,si = 1-0. The decoupling between the primary and 
the secondary channels occurs at Fu = 1.75. The value 
of F]j at decoupling is lower than for a sample containing 
only two channels since the primary channel must now 
drag twice as many secondary particles. 

In Fig. [TUTb') we plot the channel velocities versus Fjj 
for a sample with Rsi.p = 1.0 but with more particles in 
S2, Rs2,si = 1.16. For Fd < 0.37 the system is in the 



11 





1 ' 1 


1 ' 1 


- 




/ IV _ 


- 


11 f 


- 


- 




^•-^ * 






III 




I \^ 






1 , 1 





0.5 



1 

R. 



1.5 



FIG. 11: The three channel dynamic phase diagram for Fd 
vs Rs2,3i in the system from Fig. [10] with Rsi,p = 1.0. The 
locations of regions I, II, III, and IV are marked. Peaks appear 
in the value of Fd at which region I ends for the commensurate 
ratios of Rs2.si ~ 1-0 and Rs2.si ~ 2.0 
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completely locked region I, while for 0.37 < Fd < 1.625 
the particles in p and si remain locked but the particles 
in S2 partially decouple. We term this range of Fjj region 
III, and in this region Vs2 still increases with increasing 
Fd. For 1.625 < Fd < 3.56, all three of the channels are 
unlocked; we call this region IV. Within region IV, the 
velocity curves contain numerous small steps associated 
with the intermittent coupling of the particles in si and 
S2- At the low Fjj end of region IV, Vsi and Vs2 both 
decrease with increasing Fd, but for 2.7 < Fd < 3.56, 
Vs2 begins to increase with increasing Fd until Vgi and 
Ks2 join at the recoupling transition into region II. Once 
the system is in region II, both Vsi and Vs2 decrease 
monotonically with increasing Fd- 

For samples with Rsi.p = 1.0 but with increasing 
Rs2,si, the general features of the velocity force curves 
are the same as Fig. [TUT b). but the transition into region 
II is pushed to higher Fd- This is illustrated in Fig. lTUT c) 
where we plot Vp, and Vs2 versus Fd for a system 
with Rsi,p = 1.0 and Rs2,si = 1.5. Here region II does 
not appear until Fd = 12.5. Fig. [TOlc) also shows more 
clearly the increase in Vs2 just below the onset of region 
II. For samples with Rsi.p = 1.0 and Rs2,si < 1-0, only 
regions I and II occur and the velocity force curves have 
the same form as the curves illustrated in Fig. [TUT a) . 

In Fig. [TT] we map out the dynamic phase diagram 
for a three channel system with = 1-0 and varied 

Rs2.si- The value of Fd at which a transition out of re- 
gion I occurs shows commensurate peaks at Rs2,si = 1-0 
and Rs2,si — 2.0, while the region Ill-region IV tran- 
sition falls at a roughly constant value of Fd — 2.1. 
The region Il-region IV transition line shifts to slightly 
higher Fd with increasing Rs2.si- This trend continues 
for Fd values higher than those shown in Fig. [TTJ un- 
til at Rs2,si = 2.0 the II-IV transition drops to a value 



FIG. 12: Vp, Vsi, and Vs2 versus Fd for a three channel 
system with Rsi,p ~ 0.75 and varied Rs2,p. (a) At Rs2,p = 
1.0 there is a single transition from region I to region II. (b) 
At Rs2,p = 0.833, the single region I-region II transition is 
accompanied by an additional secondary maximum in Vsi and 
V32 centered at Fd ~ 1.75. 



of Fd — 7.5 (not shown in the figure). These results 
indicate that commensurability effects also occur in the 
moving phases at high Fd . 



B. Dynamics for Increased Incommensuration 

In Figs. [T^ and [T^ we plot Vp, Vsi, and Vs2 versus Fd 
for a three channel system with Rsi,p = 0.75 and var- 
ied Rs2.p- At Rs2.p = 1-0, shown in Fig. [T^ a). there is 
a single decoupling transition from region I to region II 
at Fd — 1.14. For Rs2.p < 1-0, only regions I and II 
occur; however, Vsi and Vs2 may contain additional fea- 
tures such as those shown in Fig. fT^ b) for Rs2,p — 0.833. 
Here the decoupling into region II occurs near Fd = 0.5 
which is significantly lower than the location of the I-II 
transition in the Rs2.p = 1-0 case. There is also a sec- 
ondary maximum in Vsi and Vs2 near Fd = 1.7 which is 
similar to the secondary maximum that appears in Vsi 
at incommensurate fillings in the two channel system. 
At Rs2,p — 1-25 in Fig. fOf a). the sample first transi- 
tions at Fd — 0.8 from region I to region III, where the 
particles in si and p are locked but the particles in S2 
are partially decoupled. At Fd = 0.86 the sample en- 
ters region IV where all the channels are unlocked. The 
III-IV transition is also marked by a change in sign of 
the slope of the velocity-force curves for both of the sec- 
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FIG. 13: Vp, Vsi, and Vs2 versus Fd for a three channel sys- 
tem with Rai,p ~ 0.75 and varied Rs2,p- (a) At Rs2,p ~ 1.25, 
region I is followed by a transition into region III. In region 
IV^, all the channels are unlocked, Vsi increases with increas- 
ing Fd, and Vs2 decreases with increasing Fd- In region IV^, 
Vsi decreases with increasing Fd and Vs2 increases with in- 
creasing Fd- There is a transition to region II at high Fd- 
(b) At i?s2,p = 1.583, there is a small window of region I at 
low Fd which is not highlighted on the figure. There is a 
transition directly from region III to region IV^, with region 
IV-^ absent. 



ondary channels. In Fig. I13f a) we divide region IV into 
two subregions. Just above the III-IV transition we have 
region IV"^ in which Vsi increases with increasing Fj) 
while Vs2 decreases with increasing Fq. In region IV^ 
this behavior is reversed and Vsi decreases while Vs2 in- 
creases with increasing F]j- When the particles in si and 
S2 recouple, a transition from region IV^ to region II oc- 
curs and both Vsi and K2 decrease with increasing Fjj- 
There are two distinct subregions of region IV only for 
1.0 < Rs2,p < 1-5. For Rs2,p > 1.5, only region IV^ 
appears, as shown in Fig. [T^b) for Rs2,p = 1.583 where 
the III-IV'^ transition falls at Fu = 1.3. There is a small 
window of region I that occurs at very low Fo which is 
not highlighted in the figure. 

Figure [Ml shows the Fjj versus i?s2,p phase diagram for 
the system in Fig. [TH with Rsi.p = 0.75 and d/a — 0.67. 
Here, the value of Fd at which region I ends passes 
through peaks at the commensurate values of Rs2,p 
0.25, 0.5, 0.75, 1.0, and 1.75, with a weaker peak at 
Rs2,p = 2.0. The pronounced peak at Rs2,p = 0.75 also 
corresponds to the commensurability condition Rs2.si = 
1.0. Region IV^ first appears for Rs2,p = 1.0 and van- 




FIG. 14: The dynamic phase diagram of Fd vs Rs2,p for the 
system in Fig. [12] with Rsi.p = 0.75. 
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FIG. 15: Vp, Vsi, and Va2 vs Fd for the three channel system 
with -Rsi.p = 1.25 and d/a = 0.67. (a) At Rs2.si ~ 0.6 only 
regions I and II are present, (b) At Ra2,ai = 0.86 the system 
enters region II more than once. 



ishes at Rs2,p = 1-5, which is the Rs2,si — 2.0 filling. 
The value of F]j at which the transition from region II 
— to region IV^ occurs increases with increasing Rs2.p, ex- 
cept at Rs2,p = 1.5 where the II-IV^ transition suddenly 
drops to a lower value of Fd- 

In Fig. [TSIal we plot Vp, Vsi, and Vs2 versus Fd for 
a three layer system with Rsi^p = 1.25 at Rs2.si = 0.6. 
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FIG. 16: Vy, Va\, and Vai versus Fb for the three channel 
system with i?si,p = 1.25. (a) At Rs2,si = 1.067, there is 
a transition from the locked region I to region IV^. This is 
followed by a transition to region IV^. Vs2 drops below zero 
in the region marked ND where negative drag occurs, (b) Vs2 
vs Fd from (a) showing the region of negative velocities as 
well as the existence of a local maximum at higher Fd- 



Only regions I and II are present, and there is an ad- 
ditional second broad maximum in Vsi and Vs2 cen- 
tered near Fjj — 1.5. In general, for Rsi^p — 1.25 and 
Rs2,si < 0.8 or Rs2.si > 1-0, only regions I and II ap- 
pear. For 0.86 < Rs2.si < 1-0, region II is broken into 
two sections by an intermediate transition to region IV, 
as shown in Fig. [TSlb) for Rs2,si — 0.86. The system 
passes from region I to region II, then enters region IV 
and finally returns to region II at high Fo ■ 

The dynamic phase diagrams presented here give a 
concise description of the velocity-force curves as the sys- 
tem parameters are varied. Such dynamic phase dia- 
grams have been widely used in studies of driven parti- 
cle systems such as vortices in type-II superconductors 
[40| : however, they have no connection with equilibrium 
phase diagrams obtained from systems in the thermody- 
namic limit. Having more than three phase transition 
lines meet in an equilibrium phase diagram would be 
highly unusual; however, in the nonequilibrium dynamic 
phase diagram, having more than three lines meet has no 
special implications since the lines do not represent true 
phase transition lines. Whether nonequilibrium systems 
can undergo true phase transitions that resemble equilib- 
rium phase transitions is currently a topic of active study 
and is beyond the scope of this manuscript to address. 



The appearance of multiple phases typically occurs when 
the si and S2 channels can become unlocked with each 
other. For Rsi.s2 < 1-0, incommensurations in the form 
of holes are present in one channel; however, the mobil- 
ity of the holes is less than that of the interstitials which 
arise when s2 > 1-0. 



C. Negative Drag 

For 1.0 < Rs2,si < 1-6 and fixed Rsi,p = 1.25, we show 
that a negative drag effect can occur for the particles in 
S2- During negative drag, the particles in S2 move in the 
direction opposite to the direction in which the particles 
in p are being driven. Negative drag has been observed 
in coupled ID wires where Wigner crystallization is ex- 
pected to occur [32]. In Fig. [THl a) we plot Vp, Vsi, and 
Vs2 for a three channel system with Rs2.si — 1.067. Here 
the sample is in the locked region I for Fu < 0.1. For 
0.1 < Fd < 1.4, region IV^ appears with all three chan- 
nels decoupled, Vsi increasing with increasing Fd, and 
Vs2 decreasing with increasing Fd- At Fd = 1.4 there is 
a cusp in both Vsi and K2 at the onset of region IV^. 
The cusp also marks the point at which Vs2 reaches its 
maximum negative value. In Fig. I16r a) this is labeled 
ND for the negative drag region, which extends from 
1.0 < i^D < 1.6. In Fig. \Wlh) we plot Vs2 alone ver- 
sus Fd for the system in Fig. [TBT a) showing the negative 
drag effect more clearly and also showing the presence of 
a local maximum in Vs2 at Fd = 4.5. Above this drive, 
Vs2 decreases with increasing Fd but remains positive. 

In Fig. [17] we plot the time dependent velocity v{t) of 
a single particle in each of of the three channels for the 
system in Fig. [11] at Fd = 1.36 where the particles in 
S2 undergo negative drag. The velocity of the particle 
in p is always positive and is composed of two frequen- 
cies. The velocity of the particle in si again shows two 
frequencies and drops below zero for a portion of each 
cycle; however, the overall time average of the velocity 
remains positive. The particle in S2 also experiences a 
combination of positive and negative velocities; however, 
the negative velocity portion of each cycle is greater than 
the positive velocity portion, and the particle takes a step 
backwards at the negative cusp in each cycle. It was pre- 
viously demonstrated that a system driven by two exter- 
nal ac drives can exhibit a ratchet effect in the absence 
of an asymmetric substrate (s^ - fs^ ]. In our three chan- 
nel system, when Np, Nsi, and A^s2 are all different, the 
dynamical potential produced by the particles in p and 
si acts effectively like two ac driving signals for the par- 
ticles in S2- In some cases, the interfering frequencies 
of these ac drives can create a local potential maximum 
in S2 that is moving in a direction opposite to Fd- As 
Fd is further increased, the different ac frequencies shift, 
increasing or decreasing the ratchet effect until for high 
enough Fd the coupling between S2 and the particles in 
the other channels becomes so weak that a ratchet effect 
can no longer occur. 
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FIG. 17: The time dependent velocity v of a single particle in 
each of the channels for the system in Fig. [15] at Fd = 1.36. 
Upper curve: p; middle curve: si; lower curve: S2. The 
velocity of the particle in p exhibits two frequencies and is 
always positive. The velocity of the particle in si also shows 
two frequencies and passes below zero for a portion of each 
cycle, but the time averaged velocity remains positive. The 
particle in S2 spends a larger fraction of each cycle moving 
in the negative direction, producing a negative time averaged 
velocity. 




FIG. 18: Vs2 vs Fd for the system in Fig.[T5]at Rs2,si = 1.13, 
1.2, 1.26, and 1.33, as labeled. The largest negative majcimum 
occurs for Rs2,si = 1.2. 



In Fig. [m we plot only the normalized velocities Vs2 
versus Fd for three channel samples with = 1-25 

and Rs2.si — 1-13, 1.2, 1.26, and 1.33. This shows how 
the magnitude and extent of the negative drag region 
changes with filling. The negative velocity is maximum 
for i?s2,si — 1-2 and gradually decreases with increas- 
ing Rs2,si- In Fig. [Tola), the plot of Fd versus Rs2,si is 
marked with the region ND where negative drag occurs. 
The dashed line indicates the location of the the transi- 
tions between region IV-^ and region IV^. This transi- 
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FIG. 19: (a) Fd vs _Rs2,si for the system in Fig. \TB The 
shaded region marked ND indicates where the negative drag 
for the particles in S2 occurs. Dashed line: The transition 
between regions IV* and IV^. The largest negative maxi- 
mum of Vs2 falls on this line, (b) 1142 1, the magnitude of the 
largest negative maximum in Vs2 in the negative drag region, 
vs Rs2,3i for the same system. 



tion also coincides with the maximum negative value of 
Vs2 for fixed Rs2,si- In Fig. [T9lb) we show |T42| taken at 
the IV^-IV^ transition as a function of Rsi,s2, showing 
that the overall maximum negative value of Vs2 occurs 
at Rsi,s2 — 1-2. 

In Fig. [20] we show the dynamic phase diagram of Fd 
versus Rs2.si for the three channel system with Rsi^p — 
1.25 and d/a = 0.67. There are peaks in the transition 
out of region I at Rs2,si = 0.4, 0.6 0.8, 1.2, and 1.8. 
These peaks correspond to Rs2,p = 0.5, 0.75, 1.0, 1.5, 
and 2.25, with the most prominent peak appearing at 
Rs2,p — 1-0. For Rs2,si < 0.8 the system exhibits only 
regions I and II, while for 0.8 < Rs2,si < 1-0, the transi- 
tion from region I to region II is followed by a transition 
into region IV^ at higher Fd- At even higher Fd > 7.7, 
not shown in the figure, the line marking the transition 
from region II to region IV^ changes curvature and ap- 
proaches Rs2,si = 1-0 with increasing Fd- As a result, 
for 0.8 < Rsi.s2 < 1-0 there is a high-drive transition 
from region IV^ back to region II (not shown) when the 
particles in si and S2 recouple, similar to the region IV- 
region II transition illustrated at high Fd in Fig. [TSTb) . 
At Rs2,si = 1-0, the dashed line indicates the transition 
from region IV^ to region IV"*. For Rs2,si > 1-5, the up- 
per region IV'^-region IV^ transition saturates to the line 
Fd = 0.18. Near i?s2,si = 1-0, the upper IV^-IV^ tran- 
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FIG. 20: The dynamic phase diagram of Fd vs Rs2,s1 for 
the three channel system with Rsi,p ~ 1.25 and d/a — 0.67. 
The prominent commensurate peak in the region I-region II 
transition at Rs2,si = 0.8 also corresponds to the commensu- 
rability condition of Rs2,p = 1.0. The dashed line indicates 
that at Rs2,si = 1.0, the system crosses from region IV^ to 
region IV'*. At higher Fd (not shown), the line marking the 
end of region IV"* approaches Rs2,si from above, and once 
it reaches J?32,si, region IV'* disappears. Also at higher Fd 
(not shown), the line marking the beginning of region IV^ 
approaches R32.31 from below, producing a transition from 
region IV^ to region II with increasing Fd- 



sition line approaches Rs2,si from above with increasing 
Fo, and when the transition reaches i?s2,si below F]j — 2 
(not shown in the figure), region IV'^ disappears. For 
Rs2,si > 1-0, we find no recoupling transition back into 
region II within the range Fjj < 15.0. Additionally, re- 
gion 111, where the particles in p and si are locked but the 
particles in S2 are unlocked, never occurs at all. We have 
performed additional simulations for varied Rsi.p > 1.0 
other than the value Rsi^p = 1.25 shown in Fig. [20l and 
find that the same sequence of regions illustrated in the 
figure appears in each case. 



D. Unlocking of the Central Channel 

Another possible dynamic phase has the particles in p 
and S2 locked with each other while the particles in si 
are unlocked. We term this region V, and expect it to 
occur when the average interaction between the particles 
in p and S2 is greater than the interaction between the 
particles in p and si even though the distance between 
si and p is shorter than the distance between S2 and p. 
In Fig. mja) we show an example of the occurrence of 
region V in a system with Rsi,p = 1.133, Rs2,si = 1.133, 
Rs2,p = 1-0, and d/a = 1.06. In this case the p and S2 
channels are commensurate. At low Fu, the system is 
in the locked phase I. As F]j increases, the particles in 
Si decouple from the particles in S2 and p, which remain 
locked to each other. This is indicated by the region in 
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FIG. 21: Vp, Vsi, and Vs2 vs Fd for a three channel system 
with Rsi,p — Rs2,si = 1.133, Rs2,p ~ 1.0, and d/a — 0.94. 
Here we observe a transition from region I to region V, where 
the particles in p and S2 remain locked to each other but the 
particles in si are unlocked. This is followed by region IV, 
when the particles in S2 unlock from the particles in p and 
Vsi increases with increasing Fd- (b) The same data plotted 
over a larger range of Fd shows that Vsi reaches a plateau at 
Fd = 1.9 and then decreases with increasing Fd- 



which Vsi splits away from Vp and Vs2 and increases at 
a diminished rate with increasing F]j- At F]j — 0.14, 
the particles in S2 also decouple from p and the system 
enters region IV, in which Vs2 monotonically decreases 
with increasing Fo- After the particles in S2 decouple 
from the particles in p, the coupling between the particles 
in p and si is increased, as indicated by the increase in 
the slope of Vsi at the onset of region IV. Vsi continues 
to increase with increasing Fo throughout region IV and 
even rises above Vs2 for Fd > 0.2. In Fig. I^TT b) we plot 
the same data over a larger range of Fjj to show that Vs2 
reaches a maximum value near Fjj = 1.9 before turning 
over and beginning to decrease with increasing Fo - 

The results in Fig. [211 show that it is possible to achieve 
region V in certain situations, such as when the particles 
in p and S2 are commensurate. In general it is very dif- 
ficult to obtain region V behavior in our system. The 
coupling between the particles in p and those in S2 is 
relatively weak since the distance between p and S2 is 
equal to the screening length. As a result, particles in S2 
experience a weak interaction only with those particles 
in p that lie directly above their positions, and inter- 
act much more weakly still with the other particles in 
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FIG. 22: (a) The dynamic phase diagram of Fd vs d/a for a 
system with Rs\,p = 1.133, Rs2,p = 1-0, and Rs2,si = 0.883. 
Here region V occurs for d/a > 0.75. (b) The dynamic phase 
diagram of Fd vs i?s2,3i for a system with Ra2,p = 1.0 and 
d/a = 1.06. Commensurabihty peaks appear in the transition 
out of region I at Rs2,si ~ 1.0 and Rs2,si ~ 0.5. Region V 
appears on either side of the commensuration peak at -Rs2,3i ~ 
1.0. The dashed line indicates that at Rs2,si = 1-0, the system 
passes directly from region I to region II. 
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FIG. 23; The dynamic phase diagram of Fd vs Rp,32 = 
Np/Ns2 for a system in which Np is varied. Here -Rs2,si = 
1.133 and d/a^i = 0.833, where asi is the spacing of the par- 
ticles in si. All five regions appear as marked. Region IV can 
be subdivided into regions IV* and IV^ as discussed previ- 
ously, but for clarity this subdivision is omitted here. Com- 
mensuration peaks at the transition out of region I appear at 
Rp,32 = 1.0 and at Rp^s2 = 0.58825; the latter corresponds 
to the commensurabihty condition of Rsi,p ~ 2/3. The dot- 
ted line at -Rp,s2 = 0.882 corresponds to Rai^p = 1.0 where 
the system transitions directly from region I to region III. A 
transition from region III to region IV occurs for this filling 
at Fd = 1.25 (not shown). At Rp.s2 = 1.0 the transition 
from region I to region V is marked by the thick dashed line. 
For this filling, region V ends at Fd = 0.141 and is followed 
by region IV. The transition from region IV to region II at 
high i?p,s2 continues to rise to higher values of Fd as Rp,s2 
decreases over a range of Fd larger than shown in the figure. 



p. (Note that we do not cut off the interaction at the 
screening length, but continue to compute the weak in- 
teraction out to longer distances.) This suggests that 
for different screening lengths 1/k and interchannel dis- 
tances d the coupling between particles in p and particles 
in S2 could be enhanced, producing a more widespread 
occurrence of region V and leading to additional com- 
mensuration effects. The densities of the particles in the 
channels, and not merely the ratio of their numbers, also 
plays an important role in determining which dynamical 
regions will appear. For higher particle density (smaller 
a) , the couphngs between the particles in all the channels 
are reduced, as demonstrated for the two channel case in 
Fig. Hfb). Even if the effective coupling between the par- 
ticles in p and those in S2 is strengthened by altering the 
density of the particles in the channel, this coupling must 
still be stronger than the coupling between the particles 
in p and those in si in order for region V to appear. 

In order to understand where region V occurs as a func- 
tion of the coupling between the channels, in Fig. I^^a) 
we plot the dynamic phase diagram of F]j versus d/a 
for a system with fixed Rsi.p = 1.133, Rs2,p = 1-0, and 



Rsi,s2 = 0.883. For small d/a the system passes directly 
from region I to region II. At d/a — 0.5 a window of re- 
gion IV opens between regions I and II. Region V first 
appears at d/a — 0.75, and gradually disappears for in- 
creasing Rs2.si- We next consider the case of c?/a = 1.06 
and Rs2.p = 1-0 for varied Rs2,si, as shown in Fig.l^Wb). 
Here, there is a pronounced commensurabihty peak in 
the transition out of region I at Rs2.si = 1-0, where 
all the channels contain the same number of particles. 
At Rs2.si = 1.0 the system passes directly from region 
I to region II. In windows just below and just above 
Rs2,si = 1.0 we find that region V appears and is ac- 
companied by a transition to region IV with increasing 
Fd- The width of region V grows as Rs2,si = 1-0 is ap- 
proached from either side. For Rs2,si > 1-0, region V 
gradually decreases in size with increasing Rs2,si, while 
for Rs2.si < 0.75, region V vanishes completely. For 
0.5 < Rs2.si < 0.75 the system transitions from region I 
into region IV with increasing Fjj and eventually enters 
region II at high Fd (not shown). For Rs2,si < 0.5 there 
is a only a single transition from region I to region II. A 
second commensurate peak in the transition out of region 
I appears at Rs2,si = 0.5. 
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FIG. 24: A plot of Fd vs 7?p,s2 for the system in Fig. [21 
indicating the three regions in which negative drag of the 
particles in S2 occurs. 



E. Five Dynamical Phases 

In Fig. [53] we plot the dynamic phase diagram of 
versus i?p,s2 = Np/Ns2 for a system which exhibits ah five 
phases as well as several regions where a negative drag 
effect occurs. Here we vary Np and fix Rs2,si = 1.133 and 
d/ttsi = 0.833, where Osi is the spacing of the particles 
in si. For this choice of parameters, we observe region 
V only at Rp^s2 = 1-0, the value shown in Fig. [2TJ At 
Rp,s2 = 0.882, which also corresponds to Rsi^p = 1.0, 
there is a single transition from region I to region III, 
indicated by the dashed line. Here the particles in p 
and Si are locked because they are commensurate. The 
dynamics for 0.6 < Rp.s2 < 1.75 is dominated by region 
IV. For 1.75 < Rp.s2 < 1.95, a transition from region IV 
to region II occurs at higher Fd. The location of this 
transition shifts to higher values of Fd as Rp^s2 drops 
below Rp^s2 = 1-95. For Rp^s2 > 1.95, the system goes 
directly into region II for finite Fd . Region II appears for 
high Rp^s2 since as Np increases, the effectiveness of the 
coupling between the particles in p and the particles in si 
and S2 decreases. As a result, even though the particles 
in si and S2 are incommensurate, the coupling between 
the primary and secondary channels eventually becomes 
so weak that the particles in si and S2 couple with each 
other and decouple from the particles in p. At i?p,s2 = 
0.58825 there is another peak in the transition out of 
region I produced by the commensurability condition of 
Rsi,p — 2/3 at this filling. For Rp^s2 < 0.6, region I 
grows in extent and there is a window of region HI which 
separates region I at low drives and region IV at higher 
drives. For Rp^s2 < 0.15, there is a single transition from 
region I directly to region IV. 

In Fig. [211 we indicate the regions in the Fd versus 
Rp,s2 plot where negative drag of the particles in S2 oc- 
curs for the system in Fig. [531 and in Fig. [55] we show 
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FIG. 25: Representative velocity force curves Vp, Vsi, and 
Va2 vs Fd from each of the three regions where negative drive 
occurs in Fig. [H (a) Rp,s2 = 0.35. (b) Rp,s2 = 0.823. (c) 
Rp,s2 = 1.53. 



representative velocity force curves for the three differ- 
ent negative drag regions. In Fig. [53]the largest region of 
negative drag occurs for 0.52 < Rp^s2 < 0.82. There is a 
small negative drag window near i?p,s2 = 0.3. We illus- 
trate a typical velocity force curve from this window in 
Fig.l^Sla) where we plot Vp, Vsi, and Vs2 for Rp,s2 — 0.35. 
Here the negative drag occurs in region IV. There are also 
a number of slip events which appear as sharp changes 
in Vsi near Fd = 1.25. For higher Fd beyond what is 
shown in the figure, Vs2 continues to increase back above 
zero, passes through a broad peak, and then slowly de- 
creases back toward zero at high Fd- In Fig. I^ST b) we 
plot the velocity force curves at Rs2,p = 0.823 where the 
system exhibits only region IV fiow. Here the maximum 
negative value of 17,2 occurs at Fd = 0.95 in the form of 
a cusp which is accompanied by a cusplike peak in Vsi- 
For < Fd < 0.96, l^i increases linearly with increas- 
ing Fd but the particles in si are not completely locked 
with the particles in p. This corresponds to region IV"* as 
was discussed earlier; however, in the phase diagram of 
Fig. [53] we omit the distinction between regions IV"^ and 
IV^ for clarity. In Fig. we plot the velocity force 

curves at Rp^s2 — 1-53 in the third region of negative 
drag. Here we find that the magnitude of the maximum 
negative velocity in T7,2 is reduced compared to the other 
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two negative drag regions. 

The general features of the phases outhned so far also 
occur for other parameters of density and filling, indicat- 
ing that they are robust features of the system. We have 
not observed negative drag of the particles in si or p. 

V. DISCUSSION AND SUMMARY 

We investigated a simple system consisting of two or 
three coupled ID channels of particles interacting via a 
repulsive Yukawa potential where only one of the chan- 
nels is driven. For two channel systems with an equal 
number of particles in each channel, we find a single tran- 
sition from a completely locked state to a partially decou- 
pled state where particles in the secondary channel slip 
with respect to particles in the driven channel. In the 
decoupled state, the velocity of the particles in the sec- 
ondary channel gradually decreases with increasing drive 
while the velocity of the driven particles increases linearly 
with increasing drive. When the number of particles in 
the secondary channel is slightly away from commensura- 
tion with the number of particles in the primary channels, 
a two stage decoupling transition occurs where the first 
decoupling is associated with individual slips of the in- 
commensurations or vacancies in the secondary channel. 
The velocity of the particles in the secondary channel 
contimics to increase with increasing drive until the sec- 
ond decoupling transition is reached, whereupon all the 
particles in the secondary channel begin to slip. The 
driving force at which the transition from the completely 
locked to the decoupled flow occurs has peaks at integer 
commensurate ratios of the number of particles in the 
two channels as well as at certain fractional ratios such 
as 1/2 or 3/2; however, there are no peaks for low fill- 
ing ratios since the particles in the driven channel are 
effectively moving not over a fixed substrate but over 
a distortable substrate. We also observe a ratchet ef- 
fect in the two channel system where the particles in the 
secondary channel can be rectified by an asymmetric ac 
drive applied to the primary channel. This ratchet effect 
is similar to the ratchet effect found for coupled binary 
particle species where only one species is driven. 

For three channels we find that a remarkably rich vari- 
ety of dynamical phases such as coupling and decoupling 
transitions are possible and produce a variety of com- 
mensuration effects as well as pronounced signatures in 
the velocity force curves. The commensuration effects 
occur whenever the ratio of the number of particles in at 



least two of the channels is an integer or rational fraction. 
We also observe a negative drag effect for the secondary 
channel which is furthest from the driven channel. Here, 
the particles in the secondary channel move in the di- 
rection opposite to the driving direction of the primary 
channel. When the negative drag occurs, all three chan- 
nels have incommensurate fillings. The resulting multiple 
periodic forces experienced by the particle in the furthest 
secondary channel create a bi-harmonic ratchet effect of 
a type that has been observed in systems driven with 
multiple ac drives. 

Our results could be tested for colloidal particles con- 
fined to two or three channels where one of the channels 
is driven by optical means or via microfluidics. Since 
the motion of physical colloids is never perfectly one- 
dimensional, some smearing of the effects we observe 
might occur, but the general features we describe should 
be observable. A similar cixpcriment could be performed 
in a dusty plasma system with the dust particles confined 
in grooves and driven in one dimension with a laser fo- 
cused in a single plane. Some of the effects we observe 
could be relevant for certain superconducting vortex sys- 
tems in which two different types of vortices are coupled 
and one of the two vortex types is driven with an external 
current. Additionally, these effects could also be realized 
using coupled wires in which one- dimensional Wigner 
crystal states occur. The velocity-force responses that we 
predict could be a potentially powerful method for deter- 
mining whether Wigner crystals are actually present in 
the wires. It would also be interesting to study ratchet 
effects with asymmetric ac drives for three or more chan- 
nels. Here, it may be possible to induce dc currents 
flowing in different directions for different channels. Al- 
though the system we consider appears very simple, we 
have shown that it exhibits a rich variety of behaviors 
even without substrates or other complications. If a pe- 
riodic substrate were introduced in one or more of the 
channels, we expect that an even greater variety of com- 
mensuration effects and coupling between excitations in 
the channels could occur. It would also be interesting to 
consider cases where the channels are not strictly one- 
dimensional but have a finite width to allow for transi- 
tions to buckled or zig-zag states. Even for the commen- 
surate fillings, such buckling transitions could produce 
interesting new features in the drag behavior. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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